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Abstract — A technique to compute the Cumulative Distribu- 
tion Function (CDF) of the Signal-to-Interference-plus-Noise- 
Ratio (SINR) for a wireless link with a multi-antenna, Linear, 
Minimum-Mean-Square-Error (MMSE) receiver in the presence 
of interferers distributed according to a non-homogenous Poisson 
point process on the plane, and independent Rayleigh fading 
between antennas is presented. This technique is used to compute 
the CDF of the SINR for several different models of intensity 
functions, in particular, power-law intensity functions, circular- 
symmetric Gaussian intensity functions and intensity functions 
described by a polynomial in a bounded domain. Additionally it 
is shown that if the number of receiver antennas is scaled linearly 
with the intensity function, the SINR converges in probability to 
a limit determined by the "shape" of the underlying intensity 
function. This work generalizes known results for homogenous 
Poisson networks to non-homogenous Poisson networks. 

Index Terms — MMSE, Non-homogeneous, Poisson, MIMO 

I. Introduction 

Antenna arrays can improve the performance of wireless 
networks by increasing robustness through diversity, and data 
rates through spatial multiplexing, beamforming and interfer- 
ence mitigation. The performance of multi-antenna systems in 
networks depends on environmental conditions and inter-node 
distances which effect signal and interference strengths and 
thus data rates. Hence the performance of multi-antenna sys- 
tems in spatially distributed networks have received significant 
attention in the literature. 

Multi-antenna receivers in networks with uniformly random 
spatial node distribution have been studied in several works. 
Govindasamy et. al. HI used an asymptotic analysis to ap- 
proximate the spectral efficiency with MMSE receivers, and 
Jindal et. al. [2| considered a partial zero-forcing receiver 
and found that it is possible to linearly increase the area 
spectral efficiency by simultaneously increasing the number of 
antennas and density of simultaneous transmissions. Ali et. al. 
and Louie et. al. 0] found the exact CDF of the SINR with 
MMSE receivers, with the former considering single-stream 
and the latter considering multi-stream transmissions. 

In many systems however, spatial node distributions may 
not be homogenous, such as in networks with hot-spots. Inter- 
ference modeling in non-homogenous single antenna systems 
have been studied in several works such as (5), (6), Q and 

f Corresponding author. This research was supported in part by the National 
Science Foundation under Grant CCF-1 117218. 
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Fig. 1. Randomly distributed transmitters with intensity function A(r, 8) 
O.lr" 1 , with the representative receiver at the origin of the network. 



references therein. Multi-antenna systems in non-homogenous 
networks have been studied in relatively fewer works such 
as H] who considered interference-alignment in clustered 
wireless networks, [9| who considered multi-antenna systems 
in networks with Carrier-Sensing-Multiple-Access (CSMA) 
which induces correlation between actively transmitting nodes 
and [10 1 who used an asymptotic analysis to analyze the 
spectral efficiency of non-homogenous networks with linear 
MMSE receivers. 

Here, we develop a framework to characterize the SINR of 
a representative link with multiple antennas at the receiver in 
the presence of single-antenna interferers distributed according 
to a non-homogenous Poisson Point Process on the plane, 
and independent Rayleigh fading between all antennas. The 
non-homogenous node distribution is modeled by an intensity 
function A(r, 6) which controls the likelihood of nodes occur- 
ing in a small region around the point (r, 9). Since A(r, 8) 
is deterministic (which can arise if the non-homogenity is 
predictable), this model differs works such as [6 | and [8 | where 
cluster locations are random. 

This framework is used to find a closed-form expression for 
the CDF of the SINR of a link in the center of a cluster with 
a power-law distribution of node intensities (analyzed asymp- 



totically in ifTOl ) as shown in the example in Figure [T] where 
A(r, 9) = O.lr -1 . Expressions involving generalized functions 
are also provided for networks with circularly-symmetric 
Gaussian intensities and intensity functions represented by 
a polynomial within a bounded region. The latter model is 
interesting as arbitrary continuous intensity functions can be 
approximated within bounded intervals with arbitrary accuracy 
using polynomials. By fitting polynomials to experimental 
data, this result can also be used to characterize systems for 
which node distributions are not well-modeled mathematically 
but for which experimental data on node locations exist. 
Additionally, we show that if the number of receiver antennas 
L is scaled linearly with A(r, 6), the SINR converges in 
probability to a positive deterministic value which is a function 
of A(r, 0). This latter result indicates that to the extent that 
our assumptions (in particular independent Rayleigh fading) 
hold, increasing the number of antennas per receiver can help 
scale such networks. The results here are derived combining 
techniques developed for homogenous networks in [3| and 
non-homogenous network models from IflOl . 

II. System Model 

In a circular network of radius R, a receiver centered at 
the origin is surrounded by n transmitting nodes that are 
distributed independently and randomly in the circle according 
to the Probability Density Function (PDF), f r> g(r,0), which 
is related to the intensity function A(r, ff) as follows: 



/r.eM) = -A(r,e)l {0<r<R} 
H 



(1) 



where 1 



{0<r<_R} 



is the indicator function. The number of the 



transmitters n is a Poisson random variable with mean /x, 
which can be expressed as: 



R /-27T 



rA(r, 6)d6dr 



(2) 



We assume that these transmitters, also referred to as interfer- 
es, are communicating with other receivers at locations that 
do not affect our results. In addition to the interferers, a target 
transmitter is located at a fixed distance tt away from the 
receiver at the origin. 

Loss in signal power due to propagation is modeled by the 
inverse power-law model, so that the average power p from 
a node transmitting with unit power, received at a distance r 
is p = r~ a , with the path-loss exponent a > 2. The receiver 
at the origin has L antennas and each transmitter has one 
antenna. We use the subscript T to denote the target transmitter 
whereas the interferers are labeled 1,2, ■■ • , n. Using this 
notation, the transmitted symbol of the target transmitter is 
Xt and the transmitted symbol of the i-th interferer is Xi. The 
L x 1 vector gx represents the channel coefficients between 
the target transmitter and the L antennas of the receiver at the 
origin. Similarly, gi and r*j represent the channel coefficients 
and the distance between the i-th interferer and the receiver at 
the origin. The L x 1 received signal vector can be described 



by the following equation: 



-a/2 



-a/2 

gTX T + 2^ r t § 



(3) 



i=l 



where the entries of gT and gi are independent, identically- 
distributed (i.i.d.), zero-mean, unit-variance complex Gaussian 
random variables, w is a noise vector of zero-mean, i.i.d. com- 
plex Gaussian entries with variance a 2 per complex dimension. 

The receiver uses a minimum mean square error (MMSE) 
estimator to estimate xt from y. The MMSE estimator is 
known to maximize the SINR which is given by the following 
well-known formula 



SINR 



g^ (GPGt + a 2 I L ) g T 



(4) 



where the i-th column of the L x n matrix G represents the 
channel vector g^ of the i-th interferer. Il is an L x L identity 
matrix, and P = diag [r^ a , r^~ a , • • • r~ a ~\ is a diagonal matrix 
with entries corresponding to the received power from the 
interferers. Furthermore, we define the distance-normalized 
SINR as 7 — SINR • rj, to simplify notation. In deriving 
the main results we take R — > oo to model the interferers 
as resulting from a non-homogeneous Poisson point process 
with the intensity function A(r, 6). 

III. Main Results 

A. Outage Probability with General Intensity Function 

In communication systems, it is desirable to know the 
probability that the SINR is below a threshold r, which is 
often referred to as the outage probability. For a given tt, 
this probability is simply Pr{SINR < t} = F 1 (rr a ,), where 
F 7 (7) is the CDF of 7 given by the following theorem. 

Theorem 1: 



^(7) = 1 



^ W 7 )+a 2 7 ) fc 
^ fc! 

r(£>( 7 ) + a 2 7) 
T(L) 



exp(-tp(j) - cr 2 7) 



where 



A(r,fl)r T 2 a d6dr, 
1 + r Q 7 



(5) 



(6) 



and, r(.) and are the gamma function and the upper 

incomplete gamma function respectively. 
The corresponding PDF is: 



A (7) 



(^(7) + a 2 ^- 1 exgH^T) - a 2 l){<J 2 + ip'h)) 

r(L) 



where ^'(7) i s the first derivative of tpij) with respect to 7. 
Proof: Given in Appendix lAl 

As a result, increasing the number of antennas at the receiver 
from L to L + 1 reduces the outage probability by (^(7) + 
cr 2 7) i exp(-V>(7) - a 2 j)/L\. 



B. Scaling Non-homogeneous Networks by Increasing the 
Number of Antennas 

In this section, we show that the SINR on the representative 
link (link between the receiver at the origin and its target 
transmitter) converges to a constant as the density of nodes 
increases, if the number of antennas at the receiver is linearly 
increased with interferer density. This is under the assumption 
that the channel model holds (in particular the independent 
Rayleigh fading assumption) and that accurate measurements 
of channel state information are available at the receiver. A 
similar problem was addressed in the context of homogeneous 
networks in [1] and Q. 

A key result that we use to prove this is the following lemma 
which may already be known but we were not able to find it 
in the literature. 

Lemma 1: Let the upper regularized gamma function be 
denoted by Q(L, x) — r7L) , where T(L,x) is the upper 
incomplete gamma function. Let L be a positve integer and 
q > 0, then 



lim Q(L,qL) = \°> 
l^oc M if q < 1. 



(7) 



Proof: Given in Appendix IB1 

Suppose that the intensity function A(r, 9) = /3A c (r, 9), 
where A c (r, 9) is a nominal intensity function which describes 
the "shape" of the true intensity function, and (3 is the nominal 
interferer density which scales the nominal intensity function. 
We also define: 



^0(7) 



A e (r,fl)r 1 1 dffdr. (8) 
1 + r _Q 7 



Suppose that j3 = qL , i.e., the number of antennas is scaled 
linearly with the nominal density, we have ^(7) = /3t/> c (7) = 
qLip c (-y). In addition, assuming that the noise is negligible, 
the CDF in ([5]) from Theorem Q] can be expressed as: 



^(7) = 1 
By Lemma [T] we have: 



T(L,qLMl)) 
T(L) 



(9) 



lim F 7 ( 7 ) = 1 



L— »oo 



T(L,qLMl)) 

r(L) 



0, if 7 < tp~ 

1, if7>VWf 



which implies that the SINR converges in distribution to a 
constant ip' 1 f |J r^" as the number of antennas goes to 
infinity. Additionally, since convergence in distribution to a 
constant implies convergence in probability (e.g. see ifTTI ). 7 
converges in probability as well. Therefore, if we increase the 
number of antennas linearly with the nominal interferer density 
in a given network, the SINR will approach a constant non- 
zero value. This fact implies that such networks can be scaled 
by linearly increasing the number of antennas per receiver with 
user density without degrading the SINR to zero, provided that 
the assumptions are satisfied. 



C. Outage Probability with Polynomial Intensity Function 
Consider an intensity function of the following form: 



Am(r,fl) = V, a fc r fe l{o<r<-Ro} +Por e l{ % < r }, 



fe=0 



(10) 



where ao, ffli, a 2 , ...a TO are arbitrary polynomial coefficients. 
Thus, the intensity function in the range < r < Ro is 
described by a polynomial. The second term in the equation 
results in a power-law decay of the interferer density for 
r > Ro. We additionally assume that — 2 < e < — 1 and 
po > 0. Then, the corresponding CDF of 7 is: 



f£W = 1- 



where 



L-l 



i=0 



(C( 7 ) + ^ 2 7r 



exp(-^( 7 )-a 2 7) (ID 



^(7) 



n 



(2 + e)7r 



a 



+ k \ ' a ' a 7 



fe=0 



where 2-Pi(', S S ■) is Gauss's hypergeometric function. 

Since any continuous function can be uniformly approxi- 
mated by a polynomial function with arbitrary accuracy in a 
bounded interval, according to the Stone- Weierstrass Theorem 
1 12 1, we use the result above to derive the following theorem. 

Theorem 2: For any CDF ^(7) corresponding to the inten- 
sity function 



A(r,0) = h{r)l {0<r<Ro} + p a r e l {Ro < r} 



(13) 



where h(r) is any continuous function of r, there exist 
coefficients ao, 01, 02, ...a m such that 



lim F£( 7 ) = F( 7 ). 



(14) 



Proof: Given in Appendix ICl 

Theorem 2 allows us to approximate with arbitrary accuracy, 
the CDF of the SINR for any intensity function that is 
continuous in r within a finite domain using a polynomial 
expression. In particular, since efficient algorithms exist for 
fitting polynomials to real data, this technique could be useful 
to analyze networks whose geometrical characteristics are 
not easily captured by mathematical models, but for which 
experimental data on node positions are available. 

D. Outage Probability with Piecewise Power-law Function 
Consider an intensity function of the following form: 

pir ei for < r < Ri 
p 2 r e2 for Ri < r < R 2 



A(r,9) 



(15) 



,PmX 



for R m -\ < r < R n 



where efc > — 2 for each k for which Rk = 0. With this 
intensity function, in the range Rk-i < r < Rk, the intensity 
function of the interferers follows a power-law distribution 



with nominal density pt, and exponent e^. The CDF of 7 in 
this case is given by (0 with 



N 27rpii?? +£l / 2 



2 + ei 



ei 2 



m— 1 



i^ -a+e * X 



iFA 1 



-2 + a-e fc -2 + 2a-e fe 



x 2^1 1, 



Q 

-2 + a ■ 



"7^ 



7 



p2 — a- 



n 



27rp m 7 
2 — a + e m 
2 + a- 



-R, 



2-a+e r 



2-Fi 1 



m — 1 



-2 + 2a-e fci 
a 

'X 

2 + 2a - e m 



a a 
For the simplest case, consider the intensity function: 

AM) = pr\ 



(16) 



(17) 



where — 2 < e < a — 2, to prevent interference from going 
to infinity as R — > 00. This intensity function can be used to 
model a network with a dense cluster of interferers, centered 
on the representative receiver and is useful to model networks 
with hot-spots. In this case, evaluating ^(7) over a finite R, 
applying one of Euler's hypergeometric transforms 1131 . and 
taking R — > 00 yields, 



^(7) 



2tt 2 p 



7 (e+2)/a csc ( __ 



a \ a 

Substituting ([T8l into (0, we have the CDF of 7 as 



(18) 



^2ir 2 p ( e + 2\ 



exp 



2tt 2 p / e 

CSC 7T- 



7 



7 (e+2)/a + £7 2 7)i>< 
(e+2)/a _ 1 ( 19 ) 



cr 2 7 



a \ a 
E. Outage Probability with Gaussian Intensity Function 

Suppose that the receiver is located in the center of a 
cluster who intensity follows a circularly-symmetric Gaussian 
function. A(r, 9) can be written in the form of the PDF of a 
Rayleigh variable multiplied with a constant p expressed as: 



A(r,e) = p^e-^ 2 



(20) 



where v controls the width of the intensity function. For 
specific integer values of a, ^(7) can be evaluated in terms of 
generalized functions such as the hypergeomtric and Meijer-G 
functions. For instance, when a — 3, 



^(7) 



2V3n 



-G 



5 2 
2 5 



V216w 6 



0,0 



0,| 

112 
' 3' 2' 3 



(21) 



where G is the Meijer G-function. These expressions involve 
special functions which can be evaluated efficiently in most 
mathematical software packages. 



IV. Numerical Results 
A. Monte-Carlo Simulations 

In this section, we summarize numerical results for the cases 
analyzed in the previous sections. Monte-Carlo simulations 
were run for various intensity functions to validate the general 
technique in Theorem Q] 

Figure |2] shows the empirical probability density function 
(PDF) from 100,000 simulations of a wireless network with a 
circularly-symmetric Gaussian intensity function, along with a 
graph of the PDF from Theorem 1 . The additional parameters 
of the simulations are given in the caption. From the graph, it 
is clear that the simulations match the theoretical predictions. 
Additionally note that even with a large number of trials, the 
simulated PDF is not smooth which suggests that a purely 
simulation based approach to estimate the PDF of the SINR 
in this case is computationally prohibitive. Hence, the CDF 
and PDF given in terms of generalized functions through ( l2~Tb 
are useful as they can be evaluated efficiently. 
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Fig. 2. Comparison between the empirical and theoretical PDF of the SINR 
for the Gausian intensity function with mean number of interferers fj, = 1000, 
V = 500 , = 20, number of receiver antennas L = 10, a = 3, cr 2 = 

10- 14 . 

Figure [3] shows the simulated and theoretical PDF for the 
power-law intensity function A(r, 0) — ~y ■ The additional 
parameters of the simulation are shown in the figure caption. 
From the graph, it is clear that the theoretical prediction of 
Theorem 1 is accurate. 

B. Uniform Versus Clustered Networks 

We can use the power law intensity function of ( fTTI ) with 
different values of the exponent e to compare the SINRs be- 
tween networks with uniform and clustered node distributions. 

Note that if we keep the value of p fixed, different values of 
e result in radically different numbers of nodes in the vicinity 
of the receiver at the origin. To make a fairer comparison 
between different values of e, we adjust p such that the mean 
number of nodes that fall in a radius R c disk centered at 
the origin of the infinite network is fixed for the values of 
e under consideration. Figure @] shows the probability that 
the SINR is less than or equal to 10 for different values of 
e. For this plot, we have assumed that p = 3142 nodes on 
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function of SINR with the power-law intensity function A(r, 9) = °^£ 3 



The parameters used are 
100. 000 Monte-Carlo trials. 



10, L = 10, a 



10" 



and 



average in a circle of radius R c = 1000, with e varying from 
— 1 to 0, and values of p selected so that /i = 3142. The 
remaining parameters are specified in the caption. Note that 
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Fig. 4. Outage probability vs. e for power-law intensity functions with 
L=X, 2, 4, 8 and 12 receiver antennas. Representative link length rx = 5, 
mean number of interferers in a radius 1000 circle is fixed at [A = 3142, 
a = 4, and cr 2 = 10~ 12 . 



the outage probability increases significantly with clustering. 
For instance, with L = 4, the outage probability is just below 
10 -3 for e = 0, and for e = —0.5, the outage probability is 
greater than 0.1. Additionally, observe that it is possible to 
significantly reduce the outage probability by increasing the 
number of antennas at the receiver since L = 12 antennas at 
the receiver with e w —0.5 has the same outage probability as 
L = 4 in a homogenous network. 

C. Scaling node density by increasing number of antennas 

Section UlI-BI shows that if the number of antennas is scaled 
linearly with the intensity function, the SINR approaches a 
deterministic, non-zero value. To verify this result, we plotted 
the CDF of the SINR for the Gaussian intensity function with 

A(r,0) = /3A c (r,0) = ^e~ r2 / 2v2 . L = 1,5,10 and 20 



are considered with /3 increasing linearly with L in Figure 
The remaining parameters used for the plot are given in 
the caption. Figure [5] shows that as the number of interferers 
increases from 1 to 20, the CDF of SINR approaches a step 
function, i.e the SINR approaches a constant non-zero value in 
distribution implying that it converges in probability as well. 
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Fig. 3. Comparison between the empirical and theoretical probability density ° 
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Fig. 5. Cumulative distribution function of SINR (dB) with number of 
antennas increasing linearly with nominal interferer density, with v = 500 , 
r T = 20,a = 3 and ct 2 = 10" 14 . 



V. Summary and Conclusions 

A technique to compute the CDF of the SINR on a link with 
multiple receiver antenans in non-homogenous Poisson field 
of interferers is presented and used to find expressions for the 
CDF of the SINR for several representative node distributions 
including power-law intensity functions for which a closed 
form expression for the CDF is found. This result can be used 
to characterize the SINR in the center of a cluster and indicates 
that while the SINR is significantly smaller in the center of a 
dense cluster than in a homogenous network, the performance 
loss can be mitigated by using a larger number of antennas 
at the cost of additional complexity. Additionally, we showed 
that if the number of receiver antennas is scaled linearly with 
the node intensity function, the SINR converges in probability 
to a positive constant, indicating that it is possible to scale 
such networks by increasing the number of antennas provided 
that the system assumptions hold. These results are useful to 
designers of wireless systems where the spatial distribution of 
nodes are predictable, a situation that often arises in practice. 

Appendix 

A. Proof of Theorem [7] 

From [3 |, the CDF of 7 can be expressed as 



L-l 



F 7 ( 7 ) = 1 - exp(-cr 2 7) 



2„,\i-k 



(a 2 7 ) 



i=0 k=Q 



k\(i - k)\ 



PI 



1 +pj 



exp —p;Ep 



PI 



1 +P7 



(22) 



where E x represents the expectation with respect to the random 
variable x. Recall that p = r~ a and the locations of the 
transmitting nodes are characterized by the PDF f r ^(r,6). 
Consequently, for our network model we have the following 
expressions for the expectation in the previous equation: 



which shows that 



PI 



lim 

R— too 



From (HJ and ©, we have 



■0(7) = 
Substituting into d22i i yields 

L-l i 

F 7 ( 7 ) = l-cxp(- ( T 2 7)^^ 

i=0 k=0 



*7 



1 



-dOdr 



X-fy 



1 



P7 



(23) 



(24) 



(a 2 7 )'- fc 
k\{i-k)\ 



V' fc (7)exp(V'(7)) 



L-l 



(V( 7 )+a 2 7) 1 



exp(-^( 7 ) ~ cr 2 7) 



(25) 



lim Pr (Y < x) 



(26) 



(0 follows from equation (6.5.13) in lITJl . and the PDF is 
found by taking the derivative of the CDF and simplifying. 

B. Proof of Lemma [7] 

Let Xq,Xi,--- ,Xl-i be a set of i.i.d. Poisson random 
variables with mean q. Define their sum and average respec- 
tively by Y — Ylk=o an d Y — ^. By the weak law of 
large numbers Y — > q as L — > 00 in probability, which implies 
that Y — > q in distribution, i.e. 

f0, ifx<q 
I 1, if x > q. 

Since the sum of independent Poisson random variables is 
another Poisson random variable, Y is a Poisson random 
variable with mean qL. Thus, 

Q(L, qL) = Pr (Y < L) = Pr (Y < l) . (27) 

where gL) is the CDF of a Poisson random variable with 
mean qL. Taking the limit as L — >• 00 and substituting (l26l 
yields (O completing the proof. 

C. Proof of Theorem [2] 

According to the Stone-Weierstrass Theorem [12 |, for every 
8 > there exists A^(r, 6) such that for all r in [0,i?], 3 an 
integer N such that m> N implies: 

|A£(r,0)-A(r,0)|<« (28) 



Let <5i be the product of (5 and the maximum value of r ^_ r -Z , 

II a 

(29) 



max 



'7 



1 



7 



1 



(a-1) 



> 



Then, for any 5% > 0, 3 an integer N such that m> N implies 

< 



r "7 



max 



1 + r~ Q 7 
r~ Q 7 
1 + r~ a j 



|A^(r,fl)-A(r, 
|A^(r,e)-A(r, 



(30) 



1 ^ r ,_2 7 A^(r, 6>) is uniformly convergent to 



T £ r= £ 7 A(r,0) on [0.J2], and [0, 00] as A£(r,0) and A(r,0) 
are equal in (i?, 00). Consequently, we can move the limit 
outside the integrals in the following expression resulting in 

/•CO /•27T 

■0(7) 



lim A^{r,9)r— — —dOdr 
m— >oo 1 + r 



lim V'm (7)- 



(31) 



Moreover, we express (fTTl i as F^(j) = g^^-yj), and 
-P(7) — 5(^(7))- Since 5(^^(7)) i s a continuous function 
of ^(7), as m — > 00, "0m (7) - ► V ; (7)> which implies 

lim F^( 7 ) = F( 7 ) 

= 1 - E W7) + g27)i exp(-^( 7 ) - a 2 7) (32) 

i=0 

which completes the proof. 
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